Dynamic transition in vortex flow in strongly disordered Josephson junction arrays 

and superconducting thin films 
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disorder [ p8|JT9| and the plastic flow of vortices have been 
studied recently 20 1. Also, there is an intrinsic interest in 
the non-linear dynamics of J JA and their non-equilibrium 
properties, see for example 21 1. In this paper we will 
show that there is a dynamic transition in driven JJA 
with strong positional disorder. 

The current flowing in the junction between two super- 
conducting islands in a JJA is modeled as the sum of the 
Josephson supercurrent and the normal current Il9-Bjl|: 



J M (n)=/°(n)si n M (n) 



$o afl„(n) 
2itcRn dt 
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We study the dynamics of vortices in strongly disordered d = 2 Josephson junction arrays and 
superconducting films driven by a current. We find a dynamic phase transition in vortex flow at 
a current I p > I c . Below I p there is plastic flow characterized by an average- velocity correlation 
length scale in the direction of motion, which diverges when approaching I p . Above I p we find 
a moving vortex phase with homogeneous flow and short range smectic order. A finite size analysis 
shows that this phase becomes asymptotically a liquid for large length scales. 

PACS numbers: 74.60. Ge, 74.50+r, 74.60.Ec 

The study of non-equilibrium steady states of driven 
many-degrees-of-freedom systems with quenched disor- 
der are of importance in many condensed matter sys- 
tems [JI-21 . Examples of this problem are the dynamics 
of vortices in type II superconductors |QJ|] and charge 
density waves ||. For low driving forces the dynamics 
is dominated by disorder leading to a plastic flow regime 
■ On the other hand, for very large driving forces 
the randomness should be less important and all the in- 
ternal degrees of freedom will move more or less coher- 
ently as a whole HJ3]. Recently, Koshelev and Vinokur 
0j have proposed that there is an ordered moving vor- 
tex phase. However, Giamarchi and Le Doussal Q have 
shown that some modes of static disorder are still present 
in the moving system, leading to a moving Bragg glass 
(MBG) phase Q. In turn, Balents, Marchetti and Radz- 
ihovsky @] have argued that the driven state is a moving 
smectic (MS), consisting on liquid channels with trans- 
verse periodic order. Experimentally, studies of current- 
voltage characteristics [|l0 II and neutron-scattering ex- 
periments JT^ I have found a reordering of the vortex 
structure when increasing the current bias. Recently, 
Pardo et al. |i~3f| have found in decoration experiments 
that for low magnetic fields there is a MS vortex struc- 
ture while for high fields there is a MBG. Numerical sim- 
ulation studies have also found an ordering of the vortex 
system for high currents [p],p|Jl^-[r7|] . Moon et al JlJ] have 
found a MS phase in d = 2 molecular dynamics simula- 
tions with a short-range interaction potential, while in 
[jlTj a MBG phase was found in a driven XY model sim- 
ulation for d = 3 and large magnetic fields. Therefore, 
questions such us the existence and nature of a moving 
phase and which effects of the disorder remain once the 
vortices are in motion are currently under discussion. 

Most of the experimental systems mentioned above 
have the difficulty that there is no control of the nature 
and amount of disorder. Therefore the study of disor- 
dered Josephson junction arrays (JJA) becomes particu- 
larly promising here, since they can be specifically fab- 
ricated with controlled randomness ]l8| ]. For example, 
questions such as the dynamics of JJA with percolation 
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where /°(n) is the critical current of the junction between 
the sites n and n + /i in a square lattice [n = (n x , n v ), 
fj, = x, y] , -R/v is the normal state resistance and 0^ t (n) = 
0(n + fj.) - 0(n) - A M (n) = A M 0(n) - A M (n) is the gauge 
invariant phase difference with A^(n) — ^ J*n° + A-dl. 
In the presence of an external magnetic field H we have 
A M x (n) = A x (n) — A x (n + y) + A y (n + x) — A y (n) = 
2irf, f = Ha 2 /& and a is the array lattice spac- 
ing. Here we consider a distribution of critical currents 
7°(n) = Jo^(n) = Jo[l + 6 (RAN - 1/2)] with RAN 
a random uniform number in [0,1]. We take periodic 
boundary conditions (p.b.c) in both directions in the 
presence of an external current I ext in the y-direction 
in arrays with L x L junctions. The vector potential is 
taken as A M (n, t) = A^(n) — a M (i) where in the Landau 
gauge A x (n) = -2irfn y , A°(n) = and a M (t) will al- 
low for total voltage fluctuations. With this gauge the 
p.b.c. for the phases are: 9(n x + L,n y ) — 6(n x ,n y ) 
and 9(n Xl n y + L) = 9(n Xl n y ) — 2irfLn x . The condi- 
tion of a current flowing in the y- direction: X)n_^( n ) = 
I e xtL 2 5pi,y determines the dynamics of a M (t) |22]. Af- 
ter considering conservation of current, A M ■ /^(n) = 
I^(n) — / M (n — fi) = 0, we obtain the equations: 
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where 5 M (n) = <5 M (n) sin[A M 6>(n) - A°(n) - aj, we 
have normalized currents by Jo and time by tj = 
2ttcRnIo/^o, and we have defined the discrete Laplacian 
A 2 6>(n) = 6>(n+x)+6»(n-x) + 6»(n+y) + 6»(n-y)-40(n). 
These same equations represent the dynamics of a super- 
conducting thin film (STF) after discretization of a time- 
dependent London model. One starts with the current 
density as the sum of supercurrent and normal current: 

J = Js + Jw 

J = !^[^D*-(D*r*] + ^|(V0-^A) (4) 

TO* Z7TC Ot 4>o 

with D = V + if^A and *(r) = |*(r)| exp[i0(r)]. One 
takes the discretization r = (n x £, n a £) = £n with the 
rule D„tt(r) -> ±{vj>(n + /i) - exp[-if^A M (n)]*(n)}. 
After considering conservation of current, and assum- 
ing ^(n)! is quenched and depends only on disorder, 
one obtains the same equations as in (2-3). Now /^(n) 
has to be interpreted as current density normalized by 
J = 2eh\^ \ 2 /m( = <lV(87r 2 A 2 £), time normalized by 
t = c/(4™A 2 ), <V(n) = |*(n + //)||*(n)|/|*o| 2 , and the 
field density is / = H( 2 /$q = H/2ttH c2 . The T = 
dynamical equations (2-3) are solved with a second or- 
der Runge-Kutta algorithm with time step At — O.Obrj 
and integration time IOOOtj after a transient of 500tj. 
The discrete Laplacian is inverted with a fast Fourier + 
tridiagonalization algorithm as in po| . 

We consider here very strong disorder with S = 0.5. 
For STF, this corresponds to a extremely dense distri- 
bution of pinning sites with a pinning potential with a 
50% fluctuation in amplitude in the length scale of £, 
while for JJA it corresponds to a 50% fluctuation in the 
critical currents. The ground state I ex t = vortex con- 
figuration is a vortex glass with no structure in the struc- 
ture factor (no Bragg peaks). We study a magnetic field 
of / = 1/25 and system sizes of L = 50,100,150,200. 
We calculate the time average of the total voltage V = 
(v(t)} = (da y (t)/dt) (normalized by RnIq) as a function 
of I as shown in Fig. 1(a) for one sample of size L = 200. 
The error bars of V [obtained from the statistics of the 
time averaging of v(t)] are smaller than the symbol size in 
Fig. 1(a) (error in V ~ 10 -5 — 10~ 4 ). This suggests that 
the sample size is large enough to be self-averaging. The 
initial condition is a thermally quenched vortex configu- 
ration at / = 0. From this state the current / is slowly 
increased in steps of AI = 0.01 taking as initial condi- 
tion the last phase configuration of the previous current. 
We obtain similar results by slowly decreasing decreas- 
ing the current from a random vortex configuration at 
/ = 0.8. Above a critical current of I c = 0.105 (in units 
of Io) there is a non- linear onset of voltage with plas- 
tic flow of vortices. We study the time-averaged volt- 
age in the bonds parallel to the direction of the current 
drive: v y (n) = (d0 y (n,t)/dt), which is proportional to 
the average vortex speed in the direction of the Lorentz 



force. We see in Fig. 2(a) that near I c vortex flow is very 
inhomogeneous, as typical for the plastic flow regime, 
showing channels of flow for low currents. For increasing 
drives the flow becomes more homogeneous as shown in 
Fig. 2(b). We characterize the inhomogeinity of the flow 
with the correlation function for voltages along the direc- 
tion of motion: C v (x) — ■^J2 n v y( n ) v y( n + X *)~( v y) 2 - 
We see in Fig. 2(c) that the voltage correlation increases 
for increasing values of /. There is a characteristic cor- 
relation length £ v defined by C v (x) f» C v (0) exp(— #/£„). 
This voltage correlation length was proposed by Bhat- 
tacharya and Higgins jl0| as a characteristic length scale 
for the dynamics of plastic flow. We see in Fig. 1(b) that 
£„ increases with current I and it diverges (£„ 3> L) at a 
current I p (L). We have also analyzed the voltage corre- 
lation function along y, C v (y). In this case, there is al- 
ways a fast decay of C v {y) for any bias current, therefore 
above I p the voltage distribution V y (n) becomes homoge- 
neous along the x direction but it always has fluctuations 
along the y direction. This reflects the fact that for large 
drives the effect of the disorder potential becomes neg- 
ligible only along the direction of the Lorentz force, but 
along the transverse direction is still important |7|J|] . The 
size dependence of I P (L) is shown in the inset of Fig. 1(b). 
We see that in the limit of large L, I p tends to a finite 
value of I p w 0.31 in this case. This shows that in the 
thermodynamic limit there is a dynamic phase transition. 

Above I p , we study the ordering of the moving vortex 
structure. In order to follow the vortex positions directly, 
we obtain the vorticity at the plaquette n (associated to 
the site n) as 6(n) = — A p x nint[6> M (n)/27r] with nintja;] 
the nearest integer of x. We calculate the time-averaged 
vortex structure factor as 5(k) = | (-^y 5^ fi 6(n) exp(ik • 
n)| 2 . For currents I < I p , 5(k) has only the density 
peak S(k = 0) = f 2 and an isotropic ring-like structure 
as expected for plastic flow |lj-|l^]. On the other hand, 
for I > I p there are well-defined peaks in S(k) as shown 
in the surface plot of Fig. 3(a) for L = 150. We see that 
there are two strong peaks in the k y direction, at vec- 
tors K = ±Q s y, consistent with smectic ordering in the 
direction transverse to motion [H|. (Q s — 2ir/a s with 
a s w 4.5 w (v / 3/2 i f) 1 / 2 the row spacing in a triangular 
lattice). There are also small satellite peaks in the other 
K directions. The position of the peaks is better seen in 
the two dimensional gray scale plot of Fig. 3(b). We see 
that the strongest spots besides K = ±Q s y are at the 
reciprocal space K-vectors corresponding to a triangu- 
lar lattice, suggesting the presence of some orientational 
order. |23] In Fig. 3(c) we show a finite size analysis of 
S{Q S ) ~ L~ Vb . For a moving smectic with quasi-long 
range order one expects v s < 2, while v s = 2 is the value 
for a liquid. We obtain n s = 1.96 ± 0.06 w 2. There- 
fore in the thermodynamic limit vortex flow is always 
liquid-like in nature for sufficiently large length scales, 
and smectic ordering is only a short-range phenomenon 
in d = 2, at least for strong disorder. This is consistent 
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with RG calculations which find that a moving smectic 
phase might be unstable in d = 2 Q. In Fig. 1(c) we 
show the parameter of short-range smectic order S(Q S ) 
as a function of I for a lattice size L — 200. We find 
that for large drives S(Q S ) is nearly current indepen- 
dent, it decreases when decreasing / and vanishes when 
approaching I p from above. This convalidates the result 
that there is a dynamic transition at I p where there is an 
onset of anisotropic short range order. 

The "flux flow noise" has been studied in many 
current- voltage measurements ^4j,^5|. Here we ana- 
lyze the voltage noise response with the power spectrum 
P(v) = |i / Q T dtv(t) exp(i27ri<)| 2 . We find that in the 
plastic flow regime, very near to I c , the power spectrum 
shows \/v noise, see Fig. 4(a). On the other hand, for 
large drives I > I p , the power spectrum tends to a fre- 
quency independent value for v — > as shown in Fig. 4(b). 
These two types of power spectra: for / < Jp and I > I p 
have also been found in the experiments [ p4| , p5[ . 

In conclusion, we find a dynamic phase transition at a 
current I p above the critical current I c . This transition is 
between two different types of liquids which differ in their 
dynamics and the spatio-temporal nature of the flow. Be- 
low I p there is a plastic flow regime ( "turbulent" ) with an 
isotropic liquid structure. The inhomogeneity of the flow 
is characterized by a finite average-voltage correlation 
length (T^| . Above I p there is a regime with homogeneous 
flow in the direction of motion ("laminar") and with a 
structure with anisotropic short range order. Finite size 
analysis shows that moving smectic order is only a short- 
range phenomenon in d — 2 S. These moving vortex 
phases could be observable experimentally in Josephson 
junction arrays using vortex-imaging techniques. 
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FIG. 1. (a) Voltage vs. current for a 200 x 200 Josephson 
junction array (or a 200£ x 200£ superconducting thin film, 
STF) with field / = 1/25 (B/H c2 = 2tt/25 for a STF) and 
disorder S — 0.5. (b) £„: average voltage correlation length 
along the direction of the Lorentz force. Inset: finite size 
dependence of I p . (c) S(Q S ): Intensity of the Bragg peak for 
smectic order. 
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FIG. 2. Time- averaged voltage distribution v y (n) for low 
currents in the plastic flow regime. The gray scale is pro- 
portional to the voltage intensity, (a) / = 0.12/o ~ l-27 c . 
(b)J = O.I6/0 « l.GIc (c) C v (x): correlation function of the 
average voltage along the direction of vortex motion. 

FIG. 3. (a) Surface intensity plot of the structure fac- 
tor S(k) of the moving vortex system for a current 
I = 0.40 Jo ?s 4J C and L = 150. We set the central peak 
5(0) = for clarity, (b) Gray scale intensity of S(k). (c) 
Finite size analysis of the smectic order S(Q S ). 

FIG. 4. Voltage power spectrum P(u), with frequencies 
normalized by vj — 27re 'p° 7 (for a STF, vj — 4ivtj\ 2 /c). 
(a) I = 0.14I ~ 1.4J C . (b)/ = 0.40/o ~ 4.0/ c . 
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